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Combining the hydrodynamic model (Hydro code) and the transport model (PACIAE model), we 
present the Hydro-PACIAE hybrid model. We use the Hydro-PACIAE hybrid model to calculate 
Au+Au collisions at ^/snn=130 GeV and Pb+Pb collisions at -^/snn=2.76 TeV. The obtained 
pseudo-rapidity and transverse momentum distributions well reproduce the experimental data. This 
shows that the Hydro-PACIAE hybrid model is useful to describe the ultrarelativistic Heavy Ion 
Collisions at RHIC and LHC Energies. We used the hybrid model to calculate the elliptical flow for 
the Pb+Pb collisions at ^/snn=2.76 TeV. The V2 values are bigger than the experimental data at 
high pr, and will become bigger when the Hydro evolution time is longer. The results indicate that 
the selection of the transition hypersurface has significant effect on the observable quantities. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Nq, 51.20.+d 
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I. INTRODUCTION 

The nucleus-nucleus collisions at the Relativistic 
Heavy Ion Collider (RHIC) produced an initial hot and 
dense matter and this has been interpreted as a strongly 
coupled quark-gluon plasma (sQGP) [3-0] • The Large 
Hadron Collider (LHC) opens up opportunities to explore 
the quark gluon plasma (QGP) in a wider energy region. 
The hydrodynamic models are useful tools to describe 
the hot and dense matter generated in the early stage of 
ultrarelativistic heavy ion collisions. Here, interactions 
are strong and frequent, so that the microscopic mod- 
els that assume binary collisions with free propagation 
of constituents between collisions have limited validity. 
The advantage of hydrodynamic models is that one can 
vary flexibly the dissipative properties and the equation 
of state (EoS) of the matter and test its consequences on 
the reaction dynamics and the outcome. On the other 
hand, the initial and final freeze-out (FO) stages of the 
reaction are outside the domain of applicability of the 
hydrodynamic model. 

Microscopic models, e.g. the transport model PA- 
CIAE, are quite successful to describe the partonic and 
hadronic stages of the collision. Unfortunately, the trans- 
port models are usually restricted to binary collisions. 
Thus their validity is limited in the hot and dense matter 
especially when the QGP is generated. Furthermore, ex- 
plaining the QGP phase transition between the hadronic 
and the partonic phase on a microscopic level is also an 
open issue that still has to be resolved. 

Recently the so-called hybrid models, combining mi- 
croscopic and macroscopic modules, have been devel- 
oped. There are several hybrid models 0-0] • In this pa- 
per, we present a hybrid model, Hydro-PACIAE, which 
combines the hydrodynamic model (Hydro code) and the 
transport model (PACIAE model). The hybrid model 
connects the different stages of the reaction and each of 



these are matched appropriately. The initial state 
for hydrodynamic model is constructed based on the so- 
called Effective String Rope model [ill, EH , and consid- 
ering that the system is close to local equilibrium at 
high energies. The subsequent fluid dynamical model 
assumes a QGP fluid in local equilibrium. Then we con- 
sider a transition from the QGP fluid phase to the par- 
tonic/hadronic phase, which is described by a combina- 
tion of a fluid dynamical (FD) and a particle transport 
model, PACIAE. Here we try to gain more insight into 
the principles, which may govern these processes. 

We use the hybrid model, Hydro-PACIAE, first to sim- 
ulate the final hadron distribution in pseudo-rapidity and 
transverse momentum, as well as comparing the model 
results with the experimental data for Au+Au collisions 
at y / SNN=130 GeV and Pb+Pb collisions at % /snn=2.76 
TeV. Then we use the hybrid model to calculate the el- 
liptic flow in Pb+Pb collisions at v / snn=2.76 TeV. 



II. HYBRID MODEL DESCRIPTION 



The most essential observables in ultrarelativistic 
heavy ion reactions [l3j from the laboratory detector can 
help us to understand the properties of the dense and 
hot matter. In general, we expect to extract more infor- 
mation from experiments, both on the Equation of State 
(EoS) and the transport properties of matter 0, fl5| . 
We have to apply a realistic and suitable description with 
fully 3+1 dimensional dynamical evolution at all stages 
of the reaction. 
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A. hydrodynamic Model for initialization in 
Hybrid model 

The hydrodynamic stage is calculated by using the 
Particle In Cells (PIC) method which uses a mesh 
of cells with smaller, Lagrangian fluid elements in them. 
The model assumes the MIT Bag Model equation of 
state and relativistic perfect fluid. The total number of 
so called "Marker particles" (which are the Lagrangian 
cells) is about 5 million, and their number remains a con- 
stant during a collision. Marker particles are defined to 
carry the same amount of conserved baryon charge, thus 
this charge is also exactly conserved by the method. Each 
cell of the calculational mesh contains a larger number 
of marker particles, placed in the mesh randomly. The 
marker particles move among the cells of the mesh during 
the fluid dynamical development. As we know, hydrody- 
namic models assuming phase and chemical equilibrium 
during all stages of the reaction cannot account for all 
observations in general. However, they provide an effi- 
cient way to describe the nearly perfect fluid phase of 
QGP. Then, we transform the fluid in the cells into a set 
of partons or hadrons, which are then described with the 
PACIAE model. 

The quark-gluon fluid is assumed to have a small vis- 
cosity [14[ , especially in the vicinity of the critical point of 
the phase transition to/from hadronic matter 15]. This 
phase transition [i~7j has significant dynamical effect on 
the hadronization, that in local, fluid dynamical phase 
equilibrium takes a long time, which contradicts to ob- 
servations. To avoid this time delay recently rapid tran- 
sition scenarios are assumed, which can be handled well 
in the hybrid approach. 

The 3+1 dimensional, relativistic FD model, we use 
to describe energetic heavy ion reactions, is well estab- 
lished and describes the measured collective flow phe- 
nomena well [3 . In our Computational Fluid Dy- 
namics (CFD) calculation the initial state model is based 
on expanding flux tubes or streaks [111, ll2| . It is a sys- 
tem where the matter is stopped within each streak, while 
streaks expand independently of each other. Thus, this 
model is applicable streak by streak and the momentum 
of the streaks varies, especially for the peripheral streaks 
where the asymmetry between the projectile and target 
contribution to the participant matter is the biggest. So, 
the streaks at the projectile and target sides move in the 
beam (z—) direction with substantial velocity difference. 
In the initial state model [TJ E3 ^he initial transverse 
velocity is zero for all fluid elements. 

In the CFD simulations with the PIC solution method 
[IH [l9| , the equations of relativistic FD were solved for 
a perfect quark-gluon fluid. At the same time the nu- 
merical method, due to the finite grid resolution, led to 
dissipation and thus to entropy production, which was 
analyzed [20| . From the entropy production we could 
determine the corresponding "numerical viscosity", and 
this was approximately the same as the estimated, low 
viscosity of the quark gluon plasma. To avoid double 



counting, i.e. over counting, of viscous dissipation, we 
did not add additional viscous terms to our CFD model 
simulations. 

In case of heavy ion reactions the flow is not stationary 
and the shear flow geometry is only present in the initial 
state. Later, due to the large pressure of QGP the plasma 
explodes and expands almost radially in a way that the 
final flow pattern at freeze out is close to spherical, some- 
what elongated longitudinally and in the reaction plane 
(±x-direction) due to the dominant elliptic flow. 



B. Basic Description in PACIAE model 



The PACIAE model [21j-|23| is a parton and hadron 
cascade model, which is based on PYTHIA [2^ . The 
model consists of four stages: The parton initialization, 
which decomposes the fluid into partons; The post tran- 
sition partonic state consists of quarks, anti-quarks, and 
gluons; Then, in the parton evolution stage, we are 
evaluating the total and the differential cross sections 
arising from the partonic evolution, which can be sim- 
ulated by the Monte Carlo method until the parton- 
parton collisions are ceased; Finally, the partons reach 
the hadronization at the moment of freeze-out (FO). 
This is described by the phenomenological fragmenta- 
tion model (String fragmentation) and the coalescence 
model [25ll26|. After hadronization, the rescattering pro- 
duces the final hadrons, and we obtain a configuration of 
hadrons in the spatial and momentum coordinates after 
the hadronization. The detailed description for PACIAE 
is the following: 

1. PARTON INITIALIZATION: In the case if we 
start directly from two colliding nuclei, an A+A 
(A+B) collision is decomposed into binary nucleon- 
nucleon (NN) collisions based on the collision geom- 
etry, so that every collision time is calculated by as- 
suming that the nucleons propagate along straight 
line trajectories and interact with the NN inelastic 
(total) cross sections. Then the initial NN collision- 
list is constructed from these NN collisions. A NN 
collision with the earliest collision time is selected 
from the collision list, and the final state of the 
collision is obtained by the PYTHIA model with 
string fragmentation switched-off. Afterwards the 
diquarks (anti-diquarks) are broken randomly into 
quark pairs (anti-quark pairs), and one obtains a 
configuration of quarks, anti-quarks, and gluons, 
beside a few hadronic remnants for a NN collision. 

In case of the hybrid model, we generate partons 
from the fluid cells of the FD model at a transition 
hyper-surface. This is described later. 

2. PARTON EVOLUTION: The evolution is made via 
parton rescattering, with 2 — > 2, binary LO-pQCD 
differential cross sections [27|- Then the total and 
differential cross sections in the parton evolution 
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(parton rescattering) are simulated by the Monte 
Carlo method. 

3. HADRONIZATION: The parton evolution stage is 
followed by the hadronization at the moment of 
partonic freeze-out. This means no more parton 
collisions. In the PACIAE model, the Lund string 
fragmentation model and phenomenological coales- 
cence model are supplied for the hadronization of 
partons after rescattering. In case of the hybrid 
model the latter one is used. 

4. After hadronization the rescattering among pro- 
duced hadrons is dealt with the usual two- 
body collision model. Only the rescatterings 
among n, K,p, n, p(u), A, A, £, S, fi, J/^> and their 
antiparticles are considered in the calculations. An 
isospin averaged parametrization formula is used 
for the hh cross section [H, [29| . We also provide 
an option for the constant total, elastic, and inelas- 
tic cross sections (cx^if — 40 mb, — 25 mb, 
°totf = 35 mb, crJJ = 10 mb) and the assumed 
ratio of inelastic to total cross section of 0.85. 



C. Connection of Stages in Hybrid Model 

We present the connection of the hydrodynamic model 
and PACIAE model, as the transfer of the fluid cell by 
cell to partons, satisfying the conservation laws based 
on the assumption of the thermal equilibrium before and 
after the transition. The detailed steps are as follows: 

(1) The hydrodynamic model, assumes a QGP fluid, 
with an EoS given by the MIT Bag model with mass- 
less quarks and gluons and a background field. Each 
Lagrangian fluid cell has a given baryon charge, which 
are rescaling the charge density in the cells of the mesh, 
n(x, y, z). We may deal with the transition hyper-surface 
in 3 optional ways: The first case is the simplest as- 
sumption that transition takes place isochronously, i.e. 
at t = const. ; The second case is that the transition hap- 
pens on a hyperboloid, i.e. at r = \/t 2 — r 2 = const., 
where t is the proper transition time from a space-time 
point. The third case is when the transition is described 
as a sudden change across a space-time hyper-surface, 
which hyper-surface is calculated based on a condition 
(e.g. a transition temperature), if the transition condi- 
tion is known. To be simple, in this work, we choose the 
first case, in which the direction of the flow and normal 
vectors are everywhere the same, and they are perpen- 
dicular to the transition hyper-surface. 

(2) We use the perfect relativistic fluid dynamical 
equations: 

N» u = and TZ = 0. 

where u M = (7, jv), and w = e+P, T lk = w-f 2 ViVk+PSik- 
During the transition, the above dynamical equations can 



be simplified as conservation laws: 

[T^d&v] = 0, [JVdoyJ = 0. (1) 

In our simulation, the above conservation laws have 
to be satisfied [10( and the examination of the increas- 
ing entropy should be done during the whole simulation. 
For example, baryon charge, Ni, and the energy, E^, are 
calculated according to the formula for each cell on the 
transition hyper-surface: 

Ni = V C L * rH = V C L * Nfu^, (2) 

Ei = V C L J ei = V C L J u^IfV (3) 

where the index "i" runs over the different generated par- 
ton or hadron species. V^Jj is the volume of the cell in 
the local rest frame. u M is the post transition 4-velocity, 
m m = 7(1,1?), and it M = 7(1,— v), where 7 = — v 2 . 

v is the 3-velocity of the fluid cell in the partonic or 
hadronic phase. Note that this velocity after the tran- 
sition is not the same as the flow velocity of the QGP. 
The energy momentum tensor T/ 1 includes the contri- 
bution of the particles with their kinetic energies. In the 
post transition stage, particles in partonic matter include 
quarks, anti-quarks and massless gluons. The detailes of 
this transition are described in the Appendix [V] 

(3) We utilize the Stefan Boltzmann distribution to 
generate the partons, or Jiittner distribution to generate 
the hadrons, considering the strangeness suppression and 
mass difference in the phase transition. We give the ratio 
of generated parton types in equilibrium state with con- 
sidering the quark masses via the strangeness suppression 
factor e.g^s = 0.3, according to the mass difference of 
the light quarks and strange quarks at an average tran- 
sition temperature. Of course, this ratio also affects the 
momentum distribution for different parton types. 

(4) The 4-momenta of the particles are generated ac- 
cording to the Cooper-Frye formula f30j. The baryon 
chemical potentials for the given baryon number densi- 
ties are taken into account. As we mentioned above, we 
are choosing a time step (t — Constant), and dispersing 
the partons into the fluid cells with proper distributions 
after the hyper-surface. 

These steps are done for all fluid cells of the hydro 
code, which are carrying baryon number density, pressure 
and momentum. The created partons serve as parton 
initialization for the PACIAE model. 



III. ASSESSMENT IN ULTRARELATIVISTIC 
HEAVY ION COLLISIONS 

Figures [T] and [2] show the velocity and temperature 
change in the transition from the QGP fluid to partonic 
matter as function of the energy density for different tran- 
sition time steps, 3.8 fm/c and 5.0 fm/c, in the Hydro- 
PACIAE hybrid model calculation. We can see that the 
cell by cell changes from the QGP fluid to partonic mat- 
ter are modest except at low energy densities where the 
velocity and temperature change visibly. 
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FIG. 1: Velocity change in the transition as function of the 
energy density in different transition time steps, (a) 3.8 fm/c 
and (b) 5.0 fm/c, in the Hydro-PACIAE hybrid model calcu- 
lation. The cell by cell changes from the QGP (black squares) 
to the partonic matter (red dots) are modest except at low 
energy densities where the velocity visibly decreases. 



FIG. 2: Temperature change in the transition as function 
of the energy density in different transition time steps, (a) 
3.8 fm/c and (b) 5.0 fm/c, in the hybrid model calculation. 
The cell by cell changes from the QGP (black squares) to 
the partonic matter (red dots) are modest except at low en- 
ergy densities where the temperature slightly increases and 
becomes more concentrated. 



For comparing the results of Hydro-PACIAE hybrid 
model with the results of Hydro code conveniently, we 
calculated the Au+Au collisions at ^/snn=130 GeV. 
The final charged particle pseudo-rapidity distribution is 
shown in Fig. [3J where the calculated results of Hydro- 
PACIAE model is compared with the PHOBOS data [Hj]. 
The observed centrality bins correspond to a set of im- 
pact parameters, which yield to different multiplicities 
depending on the impact parameter and the FO time. 
This is illustrated in Ref. Q1] for the case of Pb+Pb 
reactions at LHC energy, and the dependence is simi- 
lar at lower energies also. In the FD model the calcu- 
lated charged particle multiplicity, N c h, is a function of 
FO time (assuming a tpo — const. FO hyper-surface), 
and of the impact parameter, b. These multiplicities can 
matched to the measured mean charged multiplicity in 
a given experimental multiplicity bin. For example, the 
impact parameter b — 5.5 fm for Au+Au collisions cor- 
responds approximately to the experimental centrality of 
20% and b = 8.2 fm for centrality of 40%. This is similar 
to the impact parameter method in transport model [32j . 

We also calculated the transverse momentum distribu- 



tion for the Au+Au collisions at ^/snn=130 GeV and 
compared the results with the STAR data 33] in Fig. @] 
The good agreement in Figs. [3] and [4] indicate that our 
Hydro-PACIAE hybrid model can successfully reproduce 
the experiment data at RHIC energy. 

Figure [5] shows the pseudo-rapidity distribution of 
Hydro-PACIAE hybrid model for the Pb+Pb collisions 
at • v /snn=2.76 TeV together with the CMS data, which 
are taken form [34j and the pseudo-rapidity range is from 
-2.5 to 2.5. We can see that in the < 2.5 area, the hy- 
brid model results are well compared with the CMS data. 
The transverse momentum distribution for the Pb+Pb 
collisions at */snn=2.76 TeV are also calculated and the 
results are shown in Fig. [5] compared with the ALICE 
data [35| . We can see that in the low and middle pr 
regions the hybrid model results well reproduce the AL- 
ICE data, but the model results are smaller in the high 
Pt region. 

We used the Hydro-PACIAE hybrid model to cal- 
culate the elliptical flow for the Pb+Pb collisions at 
v / snn=2.76 TeV. The results are compared with the AL- 
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FIG. 3: Charged particle pseudo-rapidity distribution from 
the Hydro- PACIAE hybrid model calculations, and PHOBOS 
data [3l|] in Au+Au collisions at ^/snn=130 GeV. The exper- 
imental data are given in the indicated centrality bins (full 
symbols) and for these the corresponding mean impact pa- 
rameters (open symbols) are indicated for the hybrid model 
calculations. 



FIG. 5: Charged particle pseudo-rapidity distribution from 
the Hydro- PACIAE hybrid model calculations, and CMS data 
[34| in central collisions. The experimental data are given in 
the indicated centrality bin (full symbols) and the mean im- 
pact parameters (open symbols) are indicated for the hybrid 
model calculations. 
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FIG. 4: p T distribution for the Hydro-PACIAE hybrid 
model calculation and STAR data [33J in Au+Au collisions 
at y / SNN=130 GeV, for the same centrality bins and mean 
impact parameters as in Figure [3] 



ICE dataHH IS3] in Fig. EWe see that in the low p T re- 
gion the Hydro-PACIAE hybrid model results are smaller 
than the ALICE data, but the model results are bigger 
than the ALICE data for the high pt region. This fea- 
ture originates from the hydro model |19| , mainly because 
the Hydro code in the Hydro-PACIAE hybrid model is 



ideal fluid dynamics and has little viscosity. The FO time 
is also an important quantity for the elliptical flow. We 
compared the v-i with different FO time evolution for 40% 
centrality Pb+Pb collisions at v /snn=2.76 TeV in Fig. |U 
The vi value becomes bigger when the Hydro evolution 
time increases from 0.423 to 2.54 and 5.08 fm/c. This is 
because when the Hydro evolution time is longer more 
elliptic flow develops with time in the Hydro-PACIAE 
hybrid model. 

The directed flow, t>i(?7) for Pb+Pb collisions at 
%/ s nn=2.76 TeV was also evaluated in the hybrid model, 
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FIG. 7: v 2 distribution for the Hydro-PACIAE hybrid model 
calculation, and ALICE data in central collisions and semi 
central collisions, for the same centrality bins and mean im- 
pact parameters as in Figure [5] The 0-2% ALICE data are 
taken from [H, and 10-20%, 30-40% are taken from 
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for centrality percentage bins 40% and 60%. The ob- 
tained V\{rj) showed a random statistical fluctuation 
around zero with a t>i (^-amplitude of 0.1% for \r)\ < 4, 
which increased to 0.2% for \rj\ > 4, due to the smaller- 
charged particle multiplicity at higher pseudorapidities. 
The obtained v\{pt) was zero at low px and reached to 
1% (3%) at pr — 3(4)GeV/c respectively due to statis- 
tical fluctuations with decreasing charged particle mul- 
tiplicity. Thus, the fluctuations arising from the trans- 
port model PACIAE, dominated the smaller azimuthal 
asymmetry, Vi(rj), from the initial FD stage of the hy- 
brid model. The same situation was observed for Vi(px), 
where the FD expectation for global collective flow van- 



ishes due to symmetry reasons. 

The elliptic flow, V2{n) was also dominated by strong 
statistical fluctuations. Both for the directed and elliptic 
flow the reaction plane of the evaluation was the the- 
oretical reaction plane given by the impact parameter 
vector. The initial states, arising from the FD stage of 
the model did not include any random fluctuation. These 
results are consistent with the experimental finding that 
the central and semi-central collisions are dominated by 
random fluctuations, especially by random initial state 
fluctuations, which are not correlated with the reaction 
plane. 



IV. CONCLUSIONS 

Combining the hydrodynamic model (Hydro code) and 
the transport model (PACIAE model) , we have presented 
the Hydro-PACIAE hybrid model. We used the Hydro- 
PACIAE hybrid model to calculate Au+Au collisions at 
V s nn=130 GeV and Pb+Pb collisions at ^/snn=2.76 
TeV. The pseudo-rapidity and transverse momentum dis- 
tributions well reproduce the experimental data. This 
shows that the Hydro-PACIAE hybrid model is effective 
to describe the ultrarelativistic Heavy Ion Collisions at 
RHIC and LHC Energies. We used the hybrid model to 
calculate the elliptical flow for the Pb+Pb collisions at 
\/ s nn=2.76 TeV. The v-i values are bigger than the ex- 
perimental data in the high px region, and it shows that 
the i>2 value will become bigger when the Hydro evolution 
time is longer. This result indicates significant depen- 
dence on the choice of the transition surface between the 
two components of the hybrid model. Furthermore, the 
space-time shape of the FO surface may also influence 
the final anisotropy produced by the model. 

The fluid dynamical model needs an initial stage of 
rclativistic heavy ion collisions to start, and the so-called 
Effective String Rope model was used to generate the 
initial stage for Hydro code [U, Il2l |. We also plan to 
use the transport model, PACIAE, to generate the initial 
stage for Hydro code instead of the Effective String Rope 
model. Then we will present a full transport model with 
an intermediate hydrodynamic stage to study the rela- 
tivistic heavy ion collisions. The variability of the model 
components, will enable us a bigger flexibility to study 
the role of different model assumptions, initial state and 
transition conditions on a wide set of observables. 



V. APPENDIX 

A. The transition procedure 

Here we evaluate and present the details of the tran- 
sition from the FD stage to the parton cascade stage of 
the model in case of a t — Constant transition hypersur- 
face. All components of the energy momentum tensor of 
the QGP are described in the PIC hydrodynamic model. 
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The hydro code uses the Bag Model EoS with Np = 2. 
However, some of the fluid cells with low internal energy 
would have negative pressures due to the Bag constant, 
which are set to zero in the model. Thus, we need to 
calculate the conserved quantities in the transition from 
QGP to partonic matter both for the cases of positive 
and zero pressure. First, in the QGP phase, the energy 
density component and momentum current component, 
as well as the baryon charge current are as follows: 



1 QGP 
1 QGP 
JV QGP 
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Pi 



(zqgp + Pqgp 
{eQGP + Pqgp)1qgp v QGP 

n QGPlQGP 



QGP 
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(4) 
(5) 
(6) 



where e c f is the energy density in the QGP phase in 
the calculational frame (CF). In the general case, when 
the 3- velocity of the fluid is directed arbitrarily, vqgp = 
{vQGP-x,VQ GP - y ,VQGP-z), which is known from the hy- 
dro code, then the parton side flow will keep the velocity 
3-direction in each cell because we assumed that the tran- 
sition happens at a t = Constant hyper-surface. Then 
after transition the velocity relationship between the pre- 
and post- transition velocities can be written as: 



„QGP 



,QGP 



„QGP 



v, r = 



\vqgp\ 



jv p , v 7 



,J QGP\ 



\vqgp\ 



rVp 



Here \vqgp\ 



Vy + v 2 in the QGP phase and 



the r.h.s quantities arc known from hydro code. On the 
parton side, we have 



7p 



1 



1 — 



— , or v P 



hi -I 



and we use the EoS in the QGP phase: 

eQGP + Pqgp = 4 (Pqgp + B bag ) 

a. Positive pressure in fluid cells: When the 
fluid cells have positive pressure, based on Eqs. ((H), (0 
and using the EoS, we can obtain Pqgp and eQGP, which 
are expressed by the known e c f, and vqgp (jqgp) as: 



e cf - 4B bag j QGP 
QGP = 4V — 

41 IQGP 1 



(7) 



3e c/ + ^B bagl 2 QGP v QGP 

e Q GP = 1^2 Zi ' ® 

4 7qgp 1 

Inserting these into the Eqs. (@} and ([5]), we get the 
conserved quantities T 00 and T la in terms of e c t and 
vqgp as: 

T"00 _ „ 
1 QGP — e cf 



1 QGP 



^lQGP V QGP 



—2 —{e cS -B bag ) 

* IQGP 1 



b. Zero pressure in fluid cells: In this case, we 
use the general definition, Eqs. (@1 and ((SJ with the 
condition Pqgp — for these fluid cells, and we obtain 
the simple expressions: 



1 QGP 
1 QGP 



e cf 

e c fVQGP 



This is necessary because in these cells with zero pres- 
sure, we do not use the MIT Bag model EoS, and the 
relation between T 00 and T 10 should be consistent with 
the modification. 



B. Solving the conservation laws 

The transition from QGP to partonic phase have to 
satisfy the conservation laws, Eq. (fT]), and the solution of 
this boundary condition is done in [l0|. Here we defined 
the conserved energy momentum current, A'', 

AM _ T ^ da ^ = w u n u » dai/ _ Pg^da v , 

and the baryon current j, 

j = N^do-p = n u^da^ . (9) 

With the notation, w = e + P, the above two quantities 
are a 4- vectors and an invariant scalar, which can be con- 
sidered as crossing the hyper-surface element, moreover, 
we construct two invariant scalar equations of the con- 
served energy momentum current by (i) taking its norm, 
A M A M , and (ii) taking its projection to the normal direc- 
tion, A^dcfft, they have the expressions as: 

A^A M = w(e- P){u"dcr f ,) 2 + P 2 (da' 1 da tl ) , (10) 



A^da^ = wiu^da^f - P{da> 1 d(J tl ) 



(11) 



Eliminating the term w(w M c?cr M ) 2 from Eqs. (JTUJ) and (jlip . 
we obtain that 

A M A M = (e - P)A^da^ + eP (da^da^) . (12) 

In the same way, put Eq. ((9]) into Eq. (fTTj) . then we get 



w 



f^ = A^da^ + P(da^da fl ) 



(13) 



Thus if we know the pre transition Pqgp and da^ , we 
can get the solution from Eqs. (Ti~2"]) and (|13l) to obtain the 
post transition configuration, with the equation of state. 
For example, using the simplest case in our simulation, 
where the post transition EoS is simply Pp — Ep/3, then 
from Eq. (fT2"|) we get, 

do^dop E 2 P + 2 A^da^ E P - 3 A^A^ = , (14) 

and if we assume da^ — (1,0,0,0) and — 
7(1, vqgp, 0,0), then 

A^A, = (T°° GP ) 2 - (Tq° G p) 2 and A"dcr M = T» P) 

and then Ep can be solved from the quadratic Eq. (Q~D) 
with the help of Eqs. ^ and (0. In the same way, we 
can solve np from Eq. (0, together with Eq. (|T5|) . 
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